Non-isomorphic nucleation pathways arising from morphological transitions of liquid channels 
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Motivated by unexpected morphologies of the emerging liquid phase (channels, bulges, droplets) at 
the edge of thin, melting alkane terraces, we propose a new heterogeneous nucleation pathway. The 
competition between bulk and interfacial energies and the boundary conditions determine the growth and 
shape of the liquid phase at the edge of the solid alkane terraces. Calculations and experiments reveal a 
"pre-critical" shape transition (channel-to-bulges) of the liquid before reaching its critical volume along 
a putative shape-conserving path. Bulk liquid emerges from the new shape, and depending on the degree 
of supersaturation, the new pathway may have two, one, or zero energy barriers. The findings are broadly 
relevant for many heterogeneous nucleation processes because the novel pathway is induced by common, 
widespread surface topologies (scratches, steps, etc.). 

PACS numbers: 68.08.Bc, 47.20.Dr, 64.70.dj 



Introduction. Phase transitions are influenced by the size 
or dimensionality of the system, because of contributions 
from interfacial free energies 01 Hfl. While this has long 
been known, many details are still poorly understood. Here 
we analyze the melting process of such a confined system, 
the melting of terraced solids. Our findings show for the 
first time the relevance of morphological instabilities for 
nucleation processes, in particular how pre-critical mor- 
phological transitions can lower the nucleation barriers. 

Recently, it was shown that the melting of monolayer ter- 
races of long chain alkanes at planar silica/air interfaces J3D 
leads to moving alkane drops. The drop movement was ex- 
plained by the wetting properties of the system. However, 
the emergence of the drops before they start moving has 
never been investigated. In this paper we present new ex- 
perimental data and analyze for the first time theoretically 
the onset of the drop formation. 

The emergence of the new phase is associated with a 
nucleation process. Classical nucleation theory [;4] shows 
that for alkanes the critical (heterogeneous) nucleus size is 
of order 100 nm J5fl. This is much smaller than the ob- 
served drop sizes (microns, Fig. 1). The question arises 
what mechanism leads to the unexpectedly large incipient 
droplets. We show here that the transition to the macro- 
scopically emerging and growing new phase is critically 
coupled to a morphological transition or capillary insta- 
bility from uniform to bulged liquid channels at the ter- 
race edges. The instability provides an energetically favor- 
able pathway for the phase transition. Our theory further 
predicts that this instability can only occur when the ter- 
race edge length is sufficiently large compared to the edge 
height. This conjecture is confirmed by experiments on 
the melting of alkane multilayer island terraces of varying 
thicknesses. 

Apart from addressing the relevance of morphological 
instabilities for nucleation processes, our work shows that 
surface ledges on the nanometer scale can induce droplets 
that are orders of magnitude larger, which differs from 
previous observations of channel-to-droplet instabilities on 
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FIG. 1: (color online) (a) Moving alkane droplets appear above 
the monolayer melting temperature, (b) Histogram of the droplet 
size. The size of the droplet is equivalent to the trace width. 



surface stripes Jfjfl. Long-scale instabilities have been ob- 
served in thin layers, see e.g. 0,01. but they typically in- 
volved external forces in contrast to the systems studied 
here. 

Our system has rather general properties: (i) melting 
starts at the edge, (ii) a low contact angle of the melt with 
the surface, and (iii) the contact line is pinned at the ter- 
race edge. It describes widespread situations such as the 
melting of flat (2-dimensional) solid domains or condensa- 
tion/adsorption processes on surfaces with terraces, steps, 
etc. It is also technologically relevant IsUTlll. 

Experimental preparation We prepared planar silica sur- 
faces, which were partially covered with large (tens of mi- 
crons in lateral dimension), contingent areas or terraces of 
uniform thickness that are formed by mono- or multilayers 
of solid long chain n-alkanes C36H74. The thickness of a 
solid alkane layer is approximately identical to the molec- 
ular all-trans length (~ 4.8 nm; the molecules are oriented 
upright fl^lljll ). For more details on the substrate prepa- 
ration see 1140 and ltl5ll . 

Experimental observation on monolayer melting. If the 
temperature is gradually raised above the melting point of 
the alkane monolayer T m lflr3l . eventually alkane drops ap- 
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FIG. 2: (a) Schematic diagram of the geometry: 9 is the con- 
tact angle and a is the angle between the diagonal and horizontal 
lines, (b) Possible morphologies for a given contact angle and 
increasing liquid volume: (I-II) a uniform channel with negative 
curvature, (III) a uniform channel with positive curvature, and 
(IV) a bulged channel. In (IV), the liquid bulge is not to scale. It 
is typically much larger. 



pear at the terrace edges and start to move through the solid 
monolayer area [nil. The drops have a small contact angle, 
9 < 15°. We focus here on the early stages of the drop 
formation. We observe that the drops are always already of 
micrometer size as soon as they appear at the edge (see sup- 
porting movie). This is unexpected because the step height 
is only nanometers. The drop size is measured via atomic 
force microscopy (AFM) by cooling samples with moving 
drops very quickly, thus freezing the configuration (Fig. [TJ. 
The typical width of the orifices of the paths of the drops is 
~ 1.3/xm (histogram, Fig. Q]), about 300 times larger than 
the height of the terraces ! 

Theoretical model. Fig. [2 a) shows the geometry of the 
experimental configuration of the alkane melting from the 
terrace edge with a step of height h and length L 3> h be- 
tween the smooth and homogeneous terrace and substrate 
surfaces. The liquid wets the substrate, rim, and terrace 
surfaces partially with contact angle 9. The angle a in Fig. 
[J^a) is a measure of the width of the liquid channel. 

The competition between surface and bulk energies de- 
termines how much liquid is present at temperature T and 
which morphology will be formed. The increase in the vol- 
ume of liquid alkane corresponds to the decrease in the 
volume of the solid alkane layers (the densities of liquid 
and solid alkane are approximately equal). Due to inter- 
facial energy contributions and in agreement with exper- 
imental observations, alkane melt will first appear at the 
terrace ledge. Depending on the liquid volume V, for a 
given (low) contact angle 9, there are four distinct mor- 
phologies (Fig. Elb)). (I) A channel with negative curva- 
ture at the liquid/vapour surface and an unpinned contact 
line. (II) Like (I), but with the (left segment) contact line 
pinned at the terrace edge. (Ill) A channel with positive 
curvature, pinned at the terrace edge. (IV) A channel with 
bulges (localised drops), still pinned at the terrace edge, 
with a periodic configuration in the y-direction (bulge-to- 



bulge distance A > Ay, L = nX\ Ay is the size of the 
droplet and n is an integer). The contact line depins when 
the angle of the liquid with the terrace surface exceeds the 
equilibrium contact angle. 

The interfacial energy of the system can be written as 

Ei = T, LG A LG - S LG cos 9A SL , (1) 

S is the surface tension, A the area. The subscripts L, G 
and S denote liquid, gas, and solid. The relation 

S LG cos 9 = S SG - T, SL (2) 

defines the contact angle 9. Molecular effects are ne- 
glected, since as we shall see later, the size of the resulting 
droplet is of order of microns ifTill . 

The interfacial areas of morphologies II and III, and the 
liquid volume are calculated analytically: 
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Substituting Eq. © into £Q yields the interfacial energy. 
To first approximation, we shall assume the terrace length, 
and thus the surface energy of the vertical ledge, to be con- 
stant. For morphology IV, we use a numerical method (Sur- 
face Evolver H% ). It minimises the surface energy of a 
given volume of liquid with a prescribed geometry. Sur- 
face tensions and contact angles are inputs. 

Bulk energy contributions, AfiV, result from melting 
the solid, with Afj, = /X;(T, P) - fi s (T, P), and fi t (T, P) 
and fj, s (T, P) as chemical potentials per unit volume in the 
liquid and solid phases, respectively. We estimate yfl 

A/i — —q m {T - T m )/T m , (5) 

with q m the latent heat of melting per unit volume. Posi- 
tive and negative A/i correspond to supersaturation of the 
liquid and solid, respectively. 

In the following we shall use normalized variables E = 
E/Z LG hL, V = V/h 2 L, P L = P L h/Z LG , Afi = 
Afih/T, LG , and A = X/h. 

Coexistence of channel and bulged morphologies. Sim- 
ulations reveal that for a system as shown in Fig. |2a), be- 
sides uniform channels, other isocurvature shapes are pos- 
sible (i.e., shapes with uniform pressure within the liquid). 
These shapes are characterized by bulges of equal size Ay 
that are separated by identical channel sections with dis- 
tances A > Ay, see morphology IV in Fig. [2] 

Fig. [3]a) shows the normalized interfacial energy of 
the system Ej as a function of volume V for both, uni- 
form and bulged morphologies, respectively (9 = 10° and 
A = 1600, representing our experimental observations). 
At small volume, the system is in the uniform channel mor- 
phology. The energy minimum describes a metastable uni- 
form channel with zero curvature (zero Laplace pressure). 
With increasing volume, at the minimum, the channel cur- 
vature changes from negative to positive. 

At a volume V uc the uniform morphology becomes ab- 
solutely unstable and instantaneously transforms into a 
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FIG. 3: (a) Interfacial energy as function of volume for uni- 
form and bulged morphologies (dashed line = analytical solution; 
crosses (uniform) and triangles (bulged) = numerical results). 
Here we use A = 1600 and 6 = 10°. V uc and V hc are instability 
points for the uniform and bulged morphologies respectively, (b) 
Morphology diagram showing regions where only uniform chan- 
nels (i), bulges (iii) or both morphologies (ii) can exist. The inset 
in (b) shows the bulge/drop dimension Ay — Ay/h as a function 
of the contact angle 8. 



bulged morphology (morphology IV, triangles). This tran- 
sition at Vuc is linked to a jump in energy. There is also a 
lower instability volume V& c , where starting in the bulged 
shape with V > Vb c , upon reducing the volume, the sys- 
tem will retain its shape until it transforms into a uniform 
channel at V hc . Between V uc and V hc the system behaves 
hysteretically. 

Fig. |3jb) summarizes the results of calculations as de- 
picted in Fig. [3]a) for different values of A in a morphol- 
ogy diagram. It shows regions where only uniform chan- 
nels are possible (i), both bulges and uniform channels can 
exist (ii), and only bulges are possible (iii). The bound- 
ary between regions (ii) and (iii) also agree well with an 
analytical linear instability analysis lfl5ll . which shows that 
the lowest volume instability occurs at a = 8/2 i.e., at 
Vin = 1/2 tan (0/2) + (0/2-sin(0)/2)/4 sin 4 (0/2), the 
volume corresponding to the maximum Laplace pressure 
Pl of the bulged channel, as we shall see below. 

Quite remarkably Fig. [3jb) reveals an absolute lower 
limit for A. For a bulged channel, the quasi-cylindrical 
sleeves of the channel must have the same mean curvature 
as the quasi-spherical bulge. In order to obtain a sphere 
with the same mean curvature as a cylinder, the curvature 
radius of the sphere must be twice the curvature radius of 
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FIG. 4: Melting behaviour of solid alkane terraces with heights 
of ~ 500 nm (upper row, a-c) and ~ (lower row, d-f). The 
500 nm terrace first grows a liquid channel (a) that eventually 
destabilizes into a growing bulge (b, c). The ~ 4/im terrace also 
grows a liquid channel (d). It remains a channel (e) that eventu- 
ally converges into a drop (f) without prior bulge formation. 



the cylinder. Thus, we need a certain minimal size of the 
bulge, and a certain minimal amount of volume which we 
can extract from the quasi-cylindrical sleeves to form the 
bulge. The existence of a A min also means melting will 
proceed very differently depending on the available nor- 
malized length L = L/h. Experimentally this means, the 
melting behavior of alkane terraces should depend on their 
perimeter L respectively height h. This is indeed what we 
observed Fig. [4] (see also movies in lfl5ll ). In the first row, 
L/h ~ 200 > A min , and bulged morphology appears. On 
the other hand, L/h ~ 100 < A min in the second row 
and the system melts in a channel shape. We note that en- 
ergy minimization alone favors the formation of just a sin- 
gle bulge (i.e. A = L). However, defects and dynamical 
effects may favor the formation of several droplets. De- 
tailed descriptions of Fig. [4] are given in the supporting 
information. 

The inset in Fig. [3jb) shows the corresponding bulge 
size Ay for a fixed A = 1600 (other A give similar results) 
and various contact angles when the bulge morphologies 
first become the global minima. Ay represents the mini- 
mum moving drop dimensions because if the system forms 
bulged shapes as the temperature increases (as shown be- 
low), these bulges will grow in size and eventually become 
drops. The calculated Ay ~ 600 nm for our experimental 
conditions (h ~ 5 nm, ~ 10°) agrees quite well with the 
experimentally observed minimum drop size. 

Nucleation barrier. For the nucleation and melting path- 
way, we need the total energy of the system E = Ej + 
AjiV. Stationary shapes correspond to dE/dV = 0, with 
dEj / dV = Pl , this means 

A/2 + P L = 0. (6) 

For nucleation on a uniform substrate surface, the nucle- 
ation pathway is provided by a spherical cap that is isomor- 
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FIG. 5: (a) Laplace pressure as a function of volume for uniform 
(dashed line, analytical) and bulged (crosses, Surface Evolver) 
morphologies, (b) Nucleation energy barriers as a function of 
supersaturation. Pl,mc and Pl be aTe me Laplace pressure at the 
instability points for the uniform and bulged channels. 



phically upscaled. Now, let us consider, for the moment, 
a corresponding nucleation scenario as obtained by iso- 
morphic upscaling of the uniform channel, a scenario that 
ignores the morphological instability as illustrated in Fig. 
3(a). For the uniform channel, P L = 2 sin a sin (9 — a). 
Pl is plotted in Fig. [5ja) (dashed line). The Laplace pres- 
sure has a maximum, pL.max, at a = 9/2. This implies 
(Eq. ©) a critical temperature, T in , above which there is 
no energy barrier to complete melting. T in is given by 
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We estimate T in — T m ~ 0.2° with typical values for 
C 36 H 7 4 JH. Between < P L < pL.max, there are two 
channel volumes V with the same pressure. Hence there 
are two solutions to Eq. (O for a given supersaturation 
A/2. The solution with the lower V corresponds to a local 
minimum. The larger V corresponds to a local maximum. 
Regarding the nucleation barrier, the larger V signifies the 
peak of the energy barrier for complete melting. The en- 
ergy barrier for this classical isomorphic nucleation sce- 
nario is shown as dashed line in Fig. Hb). 

The isomorphic nucleation scenario just discussed only 
applies if the uniform channel does not undergo a morpho- 
logical transition. In particular, in the situation considered 
here, a channel with A > A min will become unstable due 
to the morphological (bulging) instability discussed in the 
preceding section before it reaches its critical nucleation 
volume as a channel. A "precritical" morphological tran- 
sition opens a new and lower energy nucleation pathway. 
The Laplace pressure of the bulged shapes is plotted in Fig. 
13 a) for A = 1600 (crossed symbols). There are three dif- 
ferent regimes of supersaturation. (i) For —A/2 < pL,bc, 
the nucleation free energy landscape has two barriers. The 
first barrier corresponds to the upper instability point, at 
which the uniform channel decays into the bulged channel. 
This barrier, shown as a dashed-dotted line in Fig. Hb), 
is not a smooth maximum but instead a sawtooth-shaped 



discontinuity. The second barrier (solid line) arises from 
the bulged channel solution. Using the same argument as 
before, the peak of the barrier corresponds to the solution 
of Eq. © with P L now the Laplace pressure of the bulged 
morphology. For small positive — A/2, the second barrier is 
much larger than the first one and dominates the nucleation 
process. As we increase —A/2 towards pL,bc, the second 
barrier of the nucleation free energy is reduced and eventu- 
ally disappears. With increasing A, the value of Plm de- 
creases. It is, however, very difficult to numerically deter- 
mine its asymptotic value, (ii) For Pz bc < —A/2 < Pl, uc , 
there is only a single barrier arising from the upper insta- 
bility point of the uniform channel. This barrier disappears 
at —A/2 = Pl.uc- It is worth noting that from our lin- 
ear instability analysis lfl5ll . we find in the limit of large A, 
P L uc ~ pL,max- (iii) Finally, the nucleation free energy 
does not exhibit any barrier for — A/2 > Pl.uc- 

Discussion and Conclusion. We explain the formation of 
liquid channels at the edges of alkane terraces, their mor- 
phological transformation into bulges and the observation 
of unexpectedly large incipient alkane drops. We present 
a general nucleation and growth scenario for the melting 
of solid terraces if the liquid melt does not completely wet 
the surfaces. Promoted by interfacial energy contributions, 
melting begins already below the bulk melting tempera- 
ture T m . The liquid forms a uniform channel whose sur- 
face curvatures changes from negative to positive when the 
temperature exceeds T m . Above T m the liquid channel is 
meta-stable because the interfacial energy suppresses com- 
plete melting. If the uniform channel grew isomorphically, 
it would reach a certain critical volume corresponding to 
a maximum of the total energy. Then the continuously 
growing bulk phase emerges. We show here for the first 
time in experiment and theory that under certain condi- 
tions (sufficiently long terrace edges) the growing channel 
transforms into a bulged morphology before it reaches its 
critical nucleation volume as a channel. This "pre-critical" 
morphological transition opens a new nucleation pathway 
for the system. Its nucleation barrier is significantly lower 
than that obtained by isomorphic upscaling of the uniform 
channel geometry. This is an important result since the 
probabilities to overcome the nucleation barrier typically 
follow exponential rules oc e ^ AB / fcT . This instability as- 
pect distinguishes our work from previous "multistep" nu- 
cleation mechanisms. For examples, nucleation pathways 
with two energy barriers but no shape instabilities have 
been reported 12111 : during the crystallization of silicon, 
it was shown that the critical (crystalline) nucleus is pre- 
ceded by the formation of high-density liquid polymorph 
1 2211 . Our investigations further reveal interesting relations 
between the bulging morphology and the terrace geometry. 

Our results are quite general and relevant beyond the 
melting of alkane terraces. The terrace edge is a paradigm 
for ubiquitous surface "roughness" topologies, natural or 
artificial (scratches). The described growth scenario of 
emerging bulk phases upon melting analogously holds for 
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many condensation, adsorption, precipitation or deposition 
processes. It suggests, for instance, that a supersaturated 
vapor preferentially forms macroscopic liquid drops at a 
long edge before this happens at a short edge of the same 
height. 
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